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function h=LS_real (R, N, M, H, Nit , beta) ; 

h=zeros (N, 1) ; 

del ta= zeros (1,2); 

Q=-beta; 

d=H/2; 

for m=l:M 

delta (1) =d; delta (2) =-d; 

d=d/2; 

for it=l:Nit 
Flag=0 ; 
for p=l :N 

[val arg] =min ( [Q(p) , -Q(p) , -R(p,p) *d] ) ; 
if arg<3 
Flag=l ; 

h (p) =h (p) +delta (arg) ; 
Q=Q+delta (arg) *R (p, : ) ; 

end 

end 

if Flag==0 break; end 

end 

end 

return 
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function h=LS_cmplx (R, N, M, H, Nit , beta) ; 

h=zeros (N, 1) ; 

delta=zeros (1, 4) ; 

Q=-beta; 

d=H/2; 

for m=l : M 

delta(l)=d; delta(2)=-d; delta (3 ) =i*d; delta (4 ) =- i*d; 
d=d/2; 

for it=l:Nit 
Flag=0 ; 
for p=l :N 

[val arg] =min ( [real (Q (p) ) , -real (Q (p) ) , . . . 

imag(Q(p) ) , -imag(Q(p) ) , -R(p,p) *d] ) ; 

if arg<5 
Flag=l ; 

h (p) =h (p) +delta (arg) ; 
Q=Q+delta (arg) *R (p, : ) ; 

end 

end 

if Flag==0 break; end 

end 

end 

return 
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function h=LS_real_beta (R, N, M, H, Nit , beta) ; 
h=zeros (N, 1) ; 
delta=zeros (1,2); 
d=H/2 ; 
for m=l:M 

delta (l)=d; delta(2)=-d; 

d=d/2; 

for it=l:Nit 
Flag=0 ; 
for p=l:N 

[val arg] =min ( [-beta (p) ,beta (p) , -R (p,p) *d] ) ; 
if arg<3 
Flag=l; 

h (p) =h (p) +delta (arg) ; 
beta=beta-delta (arg) *R (p, : ) ; 

end 

end 

if Flag==0 break; end 

end 

end 
return 
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